
	// uniform distribution
    double nextRand (const double m, const double n)
    {
        return m + rand() * (n-m)/((double)RAND_MAX);
    }
